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We study the ground state of a uniform Bose gas at zero 
temperature in the Hartree-Fock-Bogoliubov (HFB) approx- 
imation. We find a solution of the HFB equations which 
obeys the Hugenholtz-Pines theorem. This solution imposes 
a macroscopic squeezing to the condensed state and as a con- 
sequence displays large particle number fluctuations. Parti- 
cle number conservation is restored by building the appropri- 
ate U(l) invariant ground state via the superposition of the 
squeezed states. The condensed particle number distribution 
of this new ground state is calculated as well as its fluctuations 
which present a normal behavior. 
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Since the recent observation of Bose-Einstein conden- 
sation in dilute gas of trapped atoms , the theoretical 
developments of dilute Bose gas have regained a consid- 
erable interest. Among all these developments, the mean 
field Hartree-Fock-Bogoliubov (HFB) approach provides 
a framework which allows to describe the properties of 
a dilute condensate in any regime of temperature Q-^)- 
But, according to previous works |2]||, the excitation 
spectrum in the approximate HFB approach is not gap- 
less in the limit of a large particle Bose gas since it does 
not fulfill the Hugenholtz-Pines (HP) theorem. At least 
one supplementary approximation should be made in or- 
der to satisfy this theorem. Among those there are the 
" Bogoliubov" approximation valid at zero tempera- 
ture and the "Popov" approximation valid at finite tem- 
perature Q. On the other hand, both these approxima- 
tions have the disadvantage that, contrary to the HFB 
one, they cannot be deduced from a variational principle. 
As a consequence, there is no stationary ground state in 
both these approximations but a ground state which ex- 
hibits a quantum phase diffusion ||. 

In this Letter, we propose to reconsider the variational 
HFB approach by investigating another kind of ground 
state solution both satisfying the HFB equations and the 
HP theorem. Solution which differs from those discussed 
previously |2|,|| can indeed be found if we admit that the 
mean value of the two atoms amplitude (ipatpa) of the 
condensed part field operator ipo is macroscopically dif- 



ferent from the mean value of the atoms number (fy'ipo). 
This is possible if the condensed part state is macroscop- 
ically squeezed in the same way as we squeezed photons 
states in quantum optics ]lO| |. 

Another complication with the HFB approach is the 
assumption that the U(l) symmetry of the ground state 
is broken plj j. Thus, we have to face the problem of deal- 
ing with states of indefinite particle number but instead 
having a well defined phase. These states are certainly 
not physical since in any experiment the Bose gas should 
conserve its total particle number. This problem is im- 
portant since the macroscopically squeezed ground state 
we propose suffers from prohibitively large particle num- 
ber fluctuations. We solve the inconvenience by taking 
as a better ground state an adequate superposition of 
squeezed states with a fixed particle number 

The effective many body Hamiltonian of the uniform 
Bose gas is given by (h = 1): 
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where %(ct) is the annihilation (creation) operator ex- 
pressed in the k momentum space with periodic bound- 
ary conditions ($(r) = exp(ifc.r)cg/v / T^), m is the 
atomic mass, Vq = J drV(r) exp(— iq.r) are the Fourier 
components of the effective interaction pseudopotential 
V(r) = (4Tra/m)5 3 (r)(d/dr)r, and a is the scattering 
length §0. 

The HFB ansatz for a homogeneous gas is l|: 
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D* is the unitary displacement operator generating a co- 
herent state in the zero momentum mode with a mean 
number N c . It has the property: 
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is the unitary squeezing operator and creates particles 
in pairs etc 1 - . Parameters controlling the popula- 
tion of pairs are denoted while 8 is an arbitrary phase 
factor. This operator has the Bogoliubov transformation 
property: 
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Sc^S* = cosh( S £)% + sinh(s E )e 2ie c t _- ; 
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The mean particle number is given by the average value 

(N) ofN = Efe4 e fe where we denotc (^) = g( e \M e )g- 
Using the transformations properties (|J) and (|J) and the 
unitary properties DD* = 1 and SS^ = 1, we get after 
a normal ordering of the operators N = (N) — J2k 
N k = S ka N c + (cosh(2sj:) - l)/2. Similarly, we get the 
following expression for the ground state energy 0: 
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sinh(2s^)/2 and where the prime in 
the sum remind us the effect of the operator (d/dr)r in 
the pseudopotential to remove ultraviolet divergences || . 
Minimizing the grand canonical Hamiltonian (H — fiN) 
where /i is the chemical potential, we obtain the HFB 
equations |||| which can be rewritten as the closed set 
of equations: 
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where we have introduced the normal and anomalous self- 
energies respectively p|,p|Jl^]: 
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Here, n = N/V is the particle density and n c — N c /V is 
the coherent particle density. 

The long wavelength excited states are given by \B) g = 

D^S^e %e ct\0) from which we can estimate using (pfjll]) the 

longwavelength excitation spectrum e ^(9\H\ff) e ^ — E g ~ 

[(k 2 /2m - n + Sn) 2 - S 2 2 ] 1/2 . In the particular HFB 
framework, the Hugenholtz-Pines (HP) theorem states 
that n — En — E12 |^]|| which is precisely the condition 
to have a gapless theory i.e. a theory in which there is 
no gap of energy between the ground state and the long- 
wavelength excited states. In that case, the excitation 
spectrum corresponds to a phonon excitation of momen- 
tum k with a velocity v — y 4iran c /m 2 . An inspection 
of (|ll]) shows that the HP theorem requires the term in 
parentheses of ( [IT| ) to be either neglected as done in the 
Bogoliubov and Popov approximation schemes Q or ei- 
ther equal to zero. 



The essential point of this paper is that the term in 
parentheses of ([ll]) can be exactly zero if we macroscop- 
ically squeeze the condensed part state. 

Indeed, we take for the parameter a negative huge 
value in such a way that n s = N s /V = (cosh(2sg) — 
l)/(2V) is a finite quantity different from zero in the 
thermodynamic limit |0|. Since by (||) we have also the 
singular expression cosh (2sg) = 1/[1 — £ 2 2 /(£n — (J,) 2 ], 
the existence of a macroscopic parameter N s respects au- 
tomatically the HP theorem in the thermodynamic limit 
pTj . Thus, defining No = ctcg, the condensed fraction 

density uq — (No) /V is not n c but the sum no = n c + n s 
and sinh(2s g )/2 = -\n s + l/(2F)]£ 12 /(£u - fi) ~ -n s . 
The density n s is still an arbitrary parameter and must 
be ajusted to verify the HP relation in (jll]). Removing 
the zero momentum mode, we can replace the sum 

by the integral j^ys J d 3 k in Eq. (|TT|) . The integral can 
be easily carried out and gives, with the help of the HP 
relation and the Eq.(0), the result n s = 8(a 3 n c /7r) 1 / 2 n c . 

Similarly defining N e = N — Nq, the same replacement 
as above gives for the density of non-zero momentum par- 
ticles n e = (N e )/V = 8(a 3 n c /7r) 1 / 2 n c /3 which is three 
times smaller than n s . Combining these results, we can 
solve (|§-[ll]) to get closed relations between the chemical 
potential, the total particle density and the ground state 



energy: 
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These are related through the thermodynamical relation 
dEg/dN = fj,. Thus, up to the first order we recover the 
classical results derived in standard textbooks . 

The physical interpretation of the squeezing of the con- 
densed part is as follows pi]. Let us define the quadra- 
ture amplitudes V = (cg + £l)/2 and Q = (ct — cg)/(2i) of 
the zero mode operators. The average value of the two 
operators are (V) = cos 9N^ 2 and (Q) = — sin ON^ 2 
]l5| while the quantum square fluctuations around their 
averages are much smaller in absolute value (S 2 V) = 
sin 2 6W s and (5 2 Q) = cos 2 6N S {8 A = A - (A)). Rep- 
resenting these results on a phase space, we see that the 
quantum fluctuations are the largest in the angular di- 
rection and microscopic in the radial one. 

Considering the potential energy part t/(ct , eg) = 

— /ictcg + §757 (ct 2 of the grand canonical Hamilto- 
nian, the distribution of the quadrature amplitudes val- 
ues prefers to spread out in a region where the poten- 
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tial energy is minimum [Tq ] . Comparing with 0, in 
the HFB approach, V is nearly a conserved quantity for 
9 = since the fluctuations of V are microscopic and 
should be rather associated with a collective motion with 
a microscopic restoring force. The commutation relations 
[Q,'P] = i/2 imply that the fluctuations of Q are huge 
but finite. 

The only difficulty in this new approach is that the 
square quantum fluctuations of the condensed part pop- 
ulation are not normal (i.e. linear in TV ) but enormous: 
(S 2 N Q ) = 2N 2 oc N 2 . This is to be compared with the 
square quantum fluctuations of the excited states popu- 
lation which are normal jL7|: 

(S 2 N e ) = ^sinh 2 (2s f )/2 = 2(7ra 3 n c ) 1/2 N c oc N (16) 

fc#0 

On the other hand, the particle number conservation im- 
poses that these fluctuations should be equal since the 
population distributions are constrained by the relation 
N Q + N e = N. 

This apparent paradox can be easily understood if 
one remembers that \9) g is not an eigenstate of the to- 
tal number operator N but has enormous fluctuations 
dominated by (5 2 No). The ansatz is not U(l) invariant 
since the phase factor 9 is fixed. Thus, in the HFB ap- 
proach, we do not have one ground state solution but 
an infinite number of degenerate ones parametrized by 
9 £ [0, 2tt[. A improved ground state solution is the 
superposition \<p) g = J Q W ^f(ff)\6) g where the function 
f(9) is determined by minimizing the energy in the \9) g 
subspace. In this way, we find the integral equation 

IT If g( e '\(H - E o)f(9)\9)g = with the eigenvalue E 
which, using the properties g (9'\9) g = g (0\9 — 9') g and 
g(9'\H\9) g = g(0\H\6 — 8') g , becomes a convolution equa- 
tion and can be solved by making an appropriate Fourier 
transform. The new solution is an eigenvector of N and 
reads: 
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in the HFB approach the symmetry is broken and the 
number fluctuations is enlarged, ultimately the ground 
state should conserve the particle number. Had we tried 
the more elaborated expression jl7| ) as the ansatz, we 
would have obtained, after a more technical calcula- 
tion, exactly the same result (|l8|). Similarly, we can 
build the number conserving excited states |^jv) e ~ 

JiT If e~ im \0) e k which are orthogonal to @ and check 
that the excitation spectrum is unchanged. 

Second, it is interesting to see how the condensed part 
is populated. This will illustrate what are the technicali- 
ties involved to handle with ( |l7|) . To this end, we notice 
that \9) g is a generating state of \<f>ft) g which further can 
be decomposed into a product of generating states, one 



for the condensed part state \Nq) = (N( 



|0),the 



other for the non-zero momentum pairs state \N e 
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The amplitude of the coefficients r^v, PN a and qN s can 
be determined through their respective generating func- 
tions: 
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where rjv is a normalisation factor. 

Let us examine the property of this improved solution. 

First, it lowers the ground state energy. Indeed, due 
to the large fluctuations, the calculations show that the 
term ^^Nq contained in the Hamiltonian gives a larger 
contribution for \9) g than for \<f>n)g- Consequently, the 
new ground state energy reads: 
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We immediatlely check that this shift in energy does not 
affect the result (|l^). This demonstrates that, although 



An exact explicit expression for pn q can be calculated 
jlOj while an approximate one for q^ c is found using the 
saddle point method fl8| . An inspection of ( po|) shows 
that the probability distribution of having No condensed 
particles out of N in the state \4>N) g is given by P(N ) = 
\pN qN~N /r N \ 2 - More explicitely, 

P(No)= m\ —) Hn ° { (-2tanh.s 5 )V2 ) 



x exp[- 



(N - (No)) 2 
2(6 2 N e ) 



(24) 



where H n (x) are the Hermite polynomials and J\f is a 
normalisation factor and (S 2 N e ) is given in Eq. (|l6|) . The 
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population distribution is not a poissonian one but a 
squeezed state one ]To[ | whose the intensity is modulated 
by a gaussian envelope centered in (No). The particle 
number fluctuations of the condensed part is given by 
(5 2 N Q ) N = g ((f> N \5 2 N Q \(f> N ) g = (6 2 N e ) N . Again using 
the saddle point method, we calculate [ p^[ : 

(S 2 N ) N = ln[ f d9e- me Go(9)G e (6 + ia)} a=Q (25) 

£ (6 2 N e ) = 2(™ 3 n c ) 1/2 A^ c (26) 

showing that the fluctuations are indeed normal. This 
result is identical to that derived previously by Giorgini 
et al. E^] by considering the condensed part as a particle 
reservoir. 

In conclusion, we have constructed a ground state 
which is stationary in the HFB approach. The key in- 
gredient is the macroscopic squeezing of the condensed 
part. This squeezing implies prohibitively large particle 
number fluctuations which are suppressed if we build the 
appropriate U(l) invariant number conserving superpo- 
sition of the squeezed states. 

An important issue is to know how the macroscopic 
squeezing is modified at non zero temperature and how 
it changes the critical temperature fill as well as other 
quantities of experimental interest p0[ . This issue can 
be addressed by means of a variational approach P]2l| 
in the same way as for fermions in the theory of super- 
conductivity j2^]. Other topics include: the investigation 
of the validity of the non-perturbative mean field result 
( |i~8| ) according to the coupling parameter na 3 Jl9| and 
an analyse of the quantum phase diffusion, in particu- 
lar, an improved calculation of the two-time correlation 
function ($(r, t)ip(f, 0)) N §,||. These further devel- 
opments would allow to find a test for the experimental 
evidence of the macroscopic squeezing. 
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